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The gravitational redshift formula is usually derived in the 
geometric optics approximation. In this note we consider an 
exact formulation of the problem in the Schwarzschild space- 
time, with the intention to clarify under what conditions this 
redshift law is valid. It is shown that in the case of shocks the 
radial component of the Poynting vector can scale according 
to the redshift formula, under a suitable condition. If that 
condition is not satisfied, then the effect of the backscattering 
can lead to significant modifications. The obtained results 
imply that the energy fiux of the short wavelength radiation 
obeys the standard gravitational redshift formula while the 
energy flux of long waves can scale differently, with redshifts 
being dependent on the frequency. 
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I. INTRODUCTION 

The gravitational redshift formula in the Schwarzschild 
spacetime pi 




(1) 



relates the frequency co' detected by an observer located 
at R with the initial frequency a; of a photon emitted 
from an areal distance a from the gravitational center. 

Formula (|^) is derived in the approximation of geo- 
metric optics (see Q, [Q, ). The aim of this paper is 
to present an exact treatment of the problem within the 
framework of the classical wave theory; the only simpli- 
fication consists in neglecting the backreaction effect. It 
is assumed that an isolated pulse of an electromagnetic 
wave is emitted outward. Its initial support is contained 
in the annulus (a, b). Initial pulses that are characterized 
by {b — a)/a << 1 will be - for the sake of brevity - re- 
ferred to as shocks. One can define a radial component of 
the energy flux Pr by a suitable spherical projection of 
the Poynting vector. It will be shown that if the relative 
width (b — a) /a is much smaller than some power of the 
relative distance (a — 2m) /a from the horizon then the 
radial energy flux of a shock conforms to a relation like 
in Eq. (0). 

The order of this paper is following. Sec. II brings the 
main result and outlines its proof. Sec. Ill is devoted 
to the estimation of a reduced electromagnetic potential. 
In Sees IV and V we derive bounds on the radiation am- 
plitudes. Sec. VI proves the desired flux relation. Last 



section is dedicated to a short summary and discussion 
on the limitations of the redshift formula. 



II. MAIN RESULTS 

The space-time geometry is defined by the 
Schwarzschild line element. 
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where i is a time coordinate, i? is a radial coordinate that 
coincides with the areal radius and dfi^ — dO^ + sin^ Odcj)^ 
is the line element on the unit sphere, < < 27r and 
< < TT. Throughout this paper G, the Newtonian 
gravitational constant, and c, the velocity of light are 
put equal to 1. We define the Regge- Wheeler coordinate 
r* ~ R -\- 2mln(2§^ — 1) and, for the sake of concise 
notation, ?y_R = 1 ^™ 



R 



The Maxwell equations read 



0, 



(3) 



where Ff^^ — d^j^A^ — d^A^ and is the electromag- 
netic potential. It is convenient to assume Aq — Q and 
the Coulomb gauge condition V^^' = 0. It is use- 
ful to follow Wheeler Q and employ the multipole ex- 
pansion. We will consider only the dipole term and, 
more specifically, choose the potential one-form A = 
yji/2sa? 6'4'(r*, t)d(l>. A similar analysis with the same 
conlusions can be done in any multipole order. 

We seek, following |^], a reduced dipole potential 
^(r*,t) in the form = ^ + 6, where S satisfies the 
dipole equation 



2 6m/ 



Here ^'(r*, t) = dofir* - 1) + ^^^^, do = dt and / is an 
arbitrary function with support in (a, b). It is well known 
that 4* solves Maxwell equations in Minkowski spacetime 
Q. / can be uniquely determined from initial data cor- 
responding to an initially outgoing radiation. Initially 
5 = do5 = 0. 

The stress-energy tensor of the electromagnetic field is 
= (1/2)(F^^F''T - (l/4)g;^F^5i^T"5) and the time-like 
translational Killing vector is denoted as C,. 
We define the projected energy density 
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The energy Eji{t) of the electromagnetic field ^' con- 
tained in the exterior of a sphere S{R) of the radius R 
reads Ef;{t) = Jj^ drp. Let n be the unit normal to 

S{R). One finds that dtER{t) — ti^Pr; thus the energy 
flux ^iR^const (here ds = ^/rjRdt is the proper time 
interval) through S{R) equals to Pr. Pr is the surface 
integral of the normal component of the Poynting vector, 



PR{R,t)= f d5(i?)r-^CM = -4^ 



ao*a^-*. (6) 



It is convenient to introduce the reduced strength field 
amplitudes h+ and 
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h+ and 7i_ represent the outgoing radiation and ingoing 
radiation, respectively. Eq. reads, in terms of the 
amplitudes h+ and h^, 



PR{R,t) 

We define 
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a quotient of the relative width of the pulse by 5-th power 
of the relative distance rj from the event horizon. This 
quantity, with this particular power of rja, is needed in 
order to prove a sufficient condition for the validity of the 
redshift formula. It is remarkable that e depends both on 
the relative width and on the relative distance from the 
Schwarzschild horizon. We comment on this point in the 
last section. 

One can prove a number of properties that are valid for 
shocks, provided that they satisfy the condition e << 1. 
These are 

i) the amplitude h^{R) is well approximated by /i^, 
where 

h+{R,t) = —{-do + dr*)(s + dof); (10) 
VR ^ ' 

ii) the function hj^(R(i)^ t) conforms to the scaling law 
rjRh+{Rit),t) = r,ah+{RiO),t^ 0); 

iii) /i_ is negligible. 

Let us define Tr^ as a null geodesic directed outward 
from the point Rq of the initial hypersurface and let 
^Ro,iR,t) be a segment of Fi^p that connects Rq and (R, t). 



Comparing the energy fiuxes through the spheres S{R) 
(where R >> 2m) and the initial S{Rq) one obtains, 
from m) and properties i) - iii), that 
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where it is assumed that R and Rq are connected by the 
null geodesic segment fRg,{R,t)- 



III. ESTIMATING THE REDUCED 
ELECTROMAGNETIC POTENTIAL 

In this section we sketch briefly derivation of bounds 
on the electromagnetic potential. There are two bounds. 
One of them (Eqs. (|l^ and (p!3[)) applies to initial data 
(or a fixed Cauchy hypersurface), and it is essentially a 
Sobolev-type inequality. The second bound (Eq. (|2l|)) 
is proven by the application of the energy method suc- 
ceeded by a Sobolev-type argument. A detailed calcula- 
tion can be found in 

One can show, in the case of initial data of compact 
support, that on the initial hypersurface 
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where A < R < b and 



F{m,y) = y -1 + 
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1 - 2m' 
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Here m = m/a and y = b/a. Eq. (gj|) is a special case of 
a result proven in 

It is useful to proceed as in |^] and define an energy 
H{R,t) of the field 5, 



H{R,t) 
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The upper integration bound is in fact finite, since the 
original support of / was finite and therefore 6 must also 
have a finite support if t < oo. This is, however, irrele- 
vant for us. It can be shown [Bl that 



y/H{at,t) < 6m / dt 




(15) 



here the external integration countour (the dt - integral) 
coincides with null geodesic Fq. The internal (dr-) inte- 
gration is done on a fixed {t = const) Cauchy hypersur- 
face. It is useful to replace the integration parametr t 



2 



by the areal radius at and to introduce the dimesionless 
variables x — at/a,y = r/a. Provided that is done, the 
insertion of ( p^ ) into (15) leads to 



IV. BOUNDS ON THE INGOING RADIATION 

The amplitude h- of the ingoing radiation can be split 
as follows 



6- 



dx 
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The internal integration can be done explicitly, so that 



dy 
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The right hand side of ( p7[ ) is bounded from above by 
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here the integration has been extended to oo and the 
Schwarz inequality has been used. This now can be cal- 
culated exactly, but the numerical factors do not really 
matter - from our point of view the only interesting fact 
is that 



' F{m,y) 
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where C is some (time-independent) constant. Later on 
we will always use the same symbol C in order to denote 
various numerical factors that appear in the course of 
calculations. 

One can derive a bound 



R 



< 
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in order to do this, observe that 



(20) 



\j;^drdr{S/r). 



The use of the Schwarz 

inequality yields now immediately, taking into account 

(0) , ^ < Jj^dr/ir^2mrV2H - ■ /IS 



Formulae (|l|) and (1^ give an estimate 
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where 



and 



h- = h- + h, 



EE —{do+dr')-9 



h=—{do+dr*)5. 
VR 



(22) 



(23) 



(24) 



h- represents the amplitude related to ^ while h is 
the term induced by the backscatter. We shall bound 
h- within the main flow inside that region of the 
Schwarzschild spacetime that is bounded by null cones 
spanned by spherical bundles of geodesies and Ff,, 
t > 0. 

A straightforward algebra yields /i_ = ^(9o+9r*)^ 
The use of (|l|) gives 
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< j^vn^- (25) 



The function h vanishes identically at t = 0, by the 
definition of ^P. Below we derive a bound on the function 
h. 

The evolution equation (U) can be rewritten as 



{-do + dr*)(r]Rh 
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Define a null geodesic Fc that is directed inward from a 
point c of the initial hypersurface. Notice that h vanishes 
identically along the outgoing null geodesic F;,. Taking 
this into account one obtains, integrating (^) along Fc, 
c> 6, 



riRh{R, t) = dr 



drnf 



(27) 



The insertion of the estimate (|21|) on the function S as 
well of the bound (O) on /, allows one to deduce from 



\h\m < 



' F{m,y) 

RVa 



(28) 
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V. SCALING OF THE AMPLITUDE OF THE 
OUTGOING RADIATION 

First we show that in the case of shocks /i+ can be 
approximated by h^. Indeed, from the definitions of 
and h+ ((0) and (|To|), respectively) follows 



J_ 



(29) 



We can conclude, taking into account the bound (12) 
that 



I/I 



1 



\h+-h+\ = ^\^-,/F{rh;yj. 



(30) 



One can easily see, by inspection of formulae (H) and 



(12), that 



< C\/e. Therefore the difference in 



(30) becomes infinitesimally small if e << 1. 

The evolution equation (^) can be written in terms of 

h+ as 



(5o + 9,..)(?/fl/i+) 



6m/ 



(31) 



Let us recall, that the right hand side of (29) can be 



pointwise bounded, using ( p^ and (^. The integration 
of ( p^ ) along a null geodesic Tfj^, a < Rq <b yields now 



mh+{R{t),t)-moh+iRo) < c^^^^^^. (32) 

Va 

Fixing the energy Ea , one notices that by choosing e < < 
1 the right hand side of ( p^ can be made arbitrarily 
small (see the remark following Eq. (^o|)). Thus the 
product rjah^ is constant. In this case one clearly sees 
the manifestation of the redshift - the rescaling of the 
amplitude If i? >> 2m (a distant observer), then 
h+ (R) = rjRg h+ {R, t ~ 0). A similar fact was discussed 
earlier in the context of a massless scalar field theory mn . 



VI. CALCULATING THE FLUX. 

Recall the formula (||): 

TT 



PR{R,t) 



(33) 



The estimates of the preceding section allow us to write, 
in the case of shocks, the initial value of the flux 



Pi 



(i?o,0) = ^=({moh+?+OWF{m,y)) (34) 



and (if (r,<) e Tr^) 

Pr{R, t) = ^ UvRh+f + OWF{rh,y)). 



(35) 



One can show, as pointed in Sec. V, that F{rh,y) < 
Crj'^e. The product 77^/1+ is constant if e << 1, as shown 
in Sec V. Therefore one obtains, comparing ( ^ ) and (|3^), 
the required relation 



Pj,{R,t)^ j'^Pn{Ro,0). 
riR 



VII. CONCLUSIONS 



(36) 



We have shown that if a relative width (fo — a) /a of 
a shock is much smaller than some power of its relative 
distance r]a from the horizon of the Schwarzschild black 
hole, then the radial energy flux satisfies Eq. (|36|). The 
condition that e << 1 can be interpreted (appealing to 
the so-called similarity theorem of the Fourier transform 
theory ||l^) as the demand that the radiation is dom- 
inated by sufficiently high frequencies. Thus one infers 
that in the high frequency limit the backscatter is negligi- 
ble; invoking now to the geometric optics approximation, 
one can relate the energy fluxes with frequencies |^ and 
the standard redshift formula (|l|) follows. It is necessary 
to point out, that barring backscatter, all of the energy 
of an outgoing pulse would get to infinity, although its 
energy flux would have to obey (^6|). The only mecha- 
nism that can diminish the energy is the backscatter of 
the radiation off the curvature of the spacetime. 

The compactness condition e << 1 is sufficient but 
not necessary for the validity of the standard redshift 
formula. The necessary condition is probably weaker - 
(6 — a)/(a7y'^) << 1, where k is some number not big- 
ger than 5 and strictly bigger than 1. In the case of 
frequencies much lower than a critical value r]^l{b — a) 
the asymptotic energy fiux might well disobey the scaling 
law (^). That suggests that the effect of the backscat- 
tering can be discovered by the observation of discrete 
spectra. Any attempt to fit observed data to the simple 
scaling law of (|l|) would lead to redshifts depending on 
the frequency, in the case of a significant backscatter. 

In the present paper the consideration is focused only 
on the dipole term, but a similar analysis with the same 
conlusions can be done in any multipole order. The key 
part that would require a modification is the Sobolev type 
estimate of section III; the adaptation of the remaining 
parts of the procedure is straightforward. 
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